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The incomplete transfer of charge and the existence of random noise 
lead to the primary operational limitations of charge transfer devices. 
Owing to the signal dependence of the residual charge, which accumulates 
as a result of the incomplete transfer, signal detection with static detection 
levels becomes seriously impaired before the onset of significant signal 
attenuation or noise degradation. A scheme using dynamic detection levels 
is found to greatly extend the operational range of CTD's and achieves 
the minimum possible error rate for detecting uncorrelated charge packet 
sizes. By contrast, simple coding procedures are found to be ineffective 
in overcoming signal degradation due to incomplete transfer. Shannon's 
expression for maximum information transmission capacity is trans- 
formed into an expression for maximum information storage capacity. 
It is found that significantly larger storage capacities are possible with 
CTD's than have been achieved. 

PRELIMINARY REMARKS 

Proposition: Devices Function at Their Limits of Operation 

One has only to design or fabricate a device which "exceeds specs," 
and then, because of his success, receive a set of revised (and more 
demanding) specifications from the systems people to appreciate this 

1453 



1454 THE BELL SYSTEM TECHNICAL JOURNAL, NOVEMBER 1973 

very basic principle. Or, if one fixed his goals too low, some clever 
fellow will come along with a new twist which makes full use of the 
device and revolutionizes the industry. On the other hand, devices 
function at their limits of operation and no better. This lies behind 
the widely recognized importance of ascertaining fundamental opera- 
tional limitations in the early stages of device development. 1 

It is, therefore, the primary objective of this article to discuss the 
operational limitations of charge transfer devices (CTD's). At a time 
when the basic device structure to be developed is still uncertain, and 
with relatively little analysis of the charge transfer process, noise, error 
rates, storage capacities, etc. (and with even less experimental verifi- 
cation of these analyses), the results of such an article may seem very 
preliminary at best. In some respects this will be the case. However, 
sufficient progress has been made in understanding the operational 
features of CTD's, especially with respect to incomplete transfer and 
noise, that definite limitations can be placed on how large information 
storage capacity and how small (digital) error rates can be made. 

It is my intention to outline CTD operational limitations in rather 
general (even perhaps philosophical) terms and to refer to appendixes, 
existing articles, work in preparation, etc., for mathematical details. 
The results are optimistic to the extent that they indicate how much 
better one can hope to do with CTD's than is commonly envisioned. 
However, at the same time, the results are pessimistic in that it is not 
obvious how one is to achieve this optimum use. Implementing our 
dynamic detection scheme, which results in the minimum possible 
error rate for digital signals, represents, nonetheless, a major step 
towards operating CTD's near full capacity. 

I. INTRODUCTION 

According to a most significant theorem due to Shannon, 2 the maxi- 
mum information transmission capacity C T of any channel is deter- 
mined by the bandwidth and the signal/noise ratio of the channel. 
With a slight modification this theorem can be transformed into a 
theorem on the maximum information storage capacity C of any 
channel, in particular of a CTD. This C, then places an upper limit 
on the number of bits of information which can be stored in an un- 
regenerated section of a CTD. If we restrict consideration to codes in 
which the size of each charge packet is independent of preceding or 
subsequent packets, we can calculate a minimum error rate and 
specify an optimum detection scheme for digital signals. To obtain 
this result, some knowledge of the signal/noise ratio and the accumu- 
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lation of residual charge owing to the incomplete transfer of portions 
of every charge packet is necessary. 

In this paper, in order to establish the minimum possible error rate 
and maximum possible storage capacities attainable with CTD's, it 
will be necessary to first review the way incomplete charge transfer 
determines the attenuation-versus-frequency characteristic and the 
total bandwidth of a CTD. The role of incomplete charge transfer in 
signal degradation will then be discussed. Once the dependence of the 
accumulated residual charge on the preceding signal is understood, 
it is possible to devise a dynamic detection scheme for simple digital 
signals and to compare this with the basic method of absolute-ampli- 
tude or fixed-threshold detection (static levels). With any form of 
detection, noise will introduce errors. Following a review of CTD 
noise, error rates are discussed. It is found that in the presence of 
noise this dynamic detection scheme attains the theoretical minimum 
error rate. Finally, maximum information storage capacities are calcu- 
lated using the bandwidth and signal/noise ratio characteristic of a 
CTD. The methods used to obtain specific results are sufficiently 
general that they can be used, for example, to calculate error rates for 
nonoptimal detection schemes. 

II. REVIEW OF INCOMPLETE CHARGE TRANSFER, SIGNAL ATTENUATION, 
AND BANDWIDTH 

As with other electronic devices, signal attenuation and device 
bandwidth are extremely useful concepts with which to discuss the 
maximum storage capacity and minimum error rates characterizing 
the operational limitations of charge transfer devices. In a CTD, signal 
attenuation arises primarily from the incomplete transfer of charge 3-9 
and only secondarily from charge losses, for example, through thermal 
(or other) leakage currents. Let us first consider incomplete transfer, 
then signal attenuation, and finally bandwidth. 

2.1 Incomplete Transfer 

It is clear that the incomplete transfer of charge from one elemental 
cell to the next will lead to signal degradation. 7 - 9 The character of this 
degradation can be ascertained as follows. The charge Qi, t in the ith 
elemental cell at time t will be the charge in the previous cell during 
the previous transfer cycle (of period r ) diminished by the charge 
Q?_ u _ r „ left behind in the (i — l)th cell but increased by the charge 
left behind in the ith cell during the immediately preceding transfer 
Qi,t- To , and also less any charge lost (or gained) during the previous 
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storage, Qi-u-r.- Thus 

Q it t = Qi-i.t-T B - Ql-u-T, 4- Ql-r. - Q\-u-r„- (1) 

In general, this equation is very difficult to solve since the Q\ and Q[ are 
nonlinear functions of Qi. Nonlinearities, however, are common in 
electronic devices. Taking the usual approach, one makes a small- 
signal analysis in order to linearize the equations. Thus we write 
Qi.t = Qi.o + qi.t, where q il0 is the time-independent (dc) component 
(bias) and q i<t is the small (ac)-signal component. Similarly, Q b and Q l 
can be decomposed into dc and ac components. Substituting these 
into eq. (1) we can obtain two equations, one for the dc terms and one 
for the ac terms. The dc equation can give us the time-independent 
(dc) charge bias level at the output. While knowledge of this may be 
important in some applications, it does not lead to any significant 
operational limitations. Of greater importance is the solution of the 
equation for the ac terms. 

The equation for the ac terms is from eq. (1) 

qi+ut+r = Qi.t - 9?,« + Qt+i.t ~ Qi,f ( 2 ) 

As part of the linearization, one takes q\ = paq< and qi = p0q<. (a and 
can be calculated from the coefficients of a Taylor series expansion 
of Q b and Q l in terms of q:pa = dQ b /dQ and p/3 = dQ l /dQ , where 
p is the number of individual charge transfers within an elemental cell 
and Q is the charge to be transferred.) Substituting into (2) one obtains 
the basic equation for qi, t : 

qi+i.t+To = Qi.t - P a Q<.t + P*flV+i.< ~" P0Q<.* ( 3 ) 

which becomes upon taking the Fourier transform 

q i+l (o})e iuT " = qi(u>)(l - pa - p/3) + paq i+ i(u>). (4) 

Up to this point we have linearized eq. (1) and passed to the fre- 
quency domain, typical procedures in electrical engineering. We pro- 
ceed to solve eq. (4) by first calculating 
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Then we note that since g. + i(w)/gi(co) is independent of i, it follows 
that q N (a>)/q (co) = [gi+i(o;)/9i(a.)] JV , where AT is the number of ele- 
mentary cells in the shift register. Recognizing that g^(w)/g (w) is the 
transfer function of the shift register, H(u), one finds that 7 - 9 
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As discussed in Section 2.2, H(w) can in principle be used to deter- 
mine q N , t for any given sequence of input charge packets qt-N To , 

The first factor in H(co) is just the phase delay in the signal present 
even in the limit of perfect transfer (a = /3 = 0). The second factor 
contains a frequency-dependent attenuation and a further phase shift. 
To a good approximation we may write 

H(u) = A («)<!-■♦<->, (6) 

where the attenuation factor A (a>) is given by 

A(co) = e -"0 e -n<xU-co8 a>r ) (7) 

and where the phase factor 0(co) is given by 

<p(w) = Nwt + not sin wr . (8) 

(n = Np, the total number of charge transfers in the shift register 
from input to output.) With knowledge of the device transfer function 
H(ui), we can discuss the attenuation A(u) and then compute the 
device bandwidth. 

2.2 Attenuation 

The attenuation factor in eq. (7) can be interpreted as follows. The 
first factor results from charge loss. If a fraction of charge is lost 
with each charge transfer, after n transfers the fraction remaining is 
just (1 - 0)* w exp (-n0) (if n/3 2 « 1). Charge loss is clearly fre- 
quency independent. The second factor results from the incomplete 
transfer of charge. For cot tt 0, very-low-frequency components, the 
size of adjacent charge packets is approximately the same. Thus the 
charge incompletely transferred at site i is nearly compensated by the 
charge incompletely transferred at i + 1. \_-paq, • + paq i+1 tt in 
eq. (3).] Thus, apart from charge losses, g*n«g< and, hence, low- 
frequency components are expected to be attenuated very little. Equa- 
tion (7) bears this out. By contrast, if wt„ = 2vf/f„ « r (/« f /2 
where f = l/ To is the clock frequency), the attenuation is relatively 
large, exp(-2r?a). Again referring to eq. (3), /«/„/2 implies that 
qi,t and qt+i, t are ~180 degrees out of phase and q itl & — q i+ i,i. Thus 
contributions to incomplete transfer add (rather than compensate as 
for low frequencies) and, again ignoring charge loss, eq. (3) predicts 
an attenuation of (1 - 2a) n « exp(— 2na). Again, eq. (7) bears this 
out. For aiTo = tt/2 (/ = / /4) the attenuation is exp( — na), an inter- 
mediate case in which the phases of each successive packet differ by 
90 degrees. 
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One further point concerning incomplete transfer should be empha- 
sized. A charge packet which "loses" a fraction a. of its charge in each 
of n transfers might be expected to be attenuated by a factor of 
(1 - (x) n « exp(-na). However, eq. (7) for A(u) shows how sensitive 
the actual degradation of a packet is to the presence and nature of the 
other charge packets composing the signal. Thus considering one 
"isolated" charge packet can be very misleading. In Appendixes A 
and B we discuss examples of attenuation in the time domain, and in 
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Fig. 1 — Attenuation versus frequency for CTD's with 2na of 0.2, 2, 4, 5, 8. 
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Appendix C we consider A (cj) in more detail. A (u>) is plotted in Fig. 
1 as a function of co for several values of na. (/3 = 0.) 

2.3 Bandwidth 

In Fig. 1 we have plotted A («) for ^ / ^ f /2, although eq. (7) 
would seem to apply for ^ /. The reason for this is basic. According 
to Shannon's sampling theorem, 2 f /2, one-half the clock frequency 
(one-half the sampling frequency), is the maximum frequency of the 
signal which can be transmitted. Thus given the clock frequency f , 
the maximum bandwidth a CTD can have is / /2. 

Incomplete charge transfer clearly reduces the effective bandwidth 
of a CTD. This is evident from the attenuation plotted in Fig. 1. 
Normally one defines bandwidth by the size of the range of frequencies 
for which the attenuation A exceeds some fraction 8 < 1. A more 
convenient definition from the point of view of information trans- 
mission and storage capacity is that the bandwidth B be given by the 
following expression: 



B 
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(9) 



where I is a modified Bessel function. 10 In Fig. 2, B is plotted as a func- 
tion of na. A slowly varying function, B decreases as (f„/2) ■ (47rna) -i 
for na » 1 . Thus despite the rapid attenuation associated with 
na ^ 10 (e- 10 PS 0.5 X 10" 4 for / = / /4 and e~ 20 ^ 0.2 X 10~ 8 for 
/ = fo/2), the bandwidth is still approximately 0.09 X (/„/2), 9 
percent of its maximum value. The relative insensitivity of /3 to na 
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Fig. 2 — Bandwidth versus na for a CTD. 
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has an important effect on the maximum achievable CTD transmission 
and storage capacities as we shall see. 

III. THE ROLE OF INCOMPLETE CHARGE TRANSFER IN SIGNAL 
DEGRADATION 

Even though attenuation and bandwidth are frequency-domain 
(and hence analog) concepts, we shall see in Section VII that it is 
these quantities, along with the signal-power-to-noise-power ratio, 
which are needed in order to calculate the maximum information 
transmission and storage capacity for digital as well as analog signals 
in CTD's. Nonetheless, it is still helpful to discuss certain features of 
digital signals in the time domain in order to better appreciate certain 
operational limitations of CTD's. For the next few sections we shall 
ignore charge losses and consider only the much more important effects 
of incomplete charge transfer on the signal. 

Suppose a charge packet of initial size Qi (representing a digital 
"one") follows some sequence of charge packets either of initial size 
Qx or of initial size Q (representing a digital "zero"). If a is the coeffi- 
cient of incomplete transfer, then after a single transfer the original 
charge packet has been reduced by a factor of (1 — a), and, after n 
similar transfers, by a factor of (1 - a) n . In addition, the original 
packet picks up some charge left behind by the preceding packets. 
This residual charge we shall refer to as Q R , which is in general a func- 
tion of the preceding signal. Thus the size of the charge packet repre- 
senting the "one" at the output is given by 

Q = (1 - a)"Qi + Qr- (10) 

A brief analysis of eq. (10) reveals that the primary operational 
limitation imposed by incomplete transfer is the dependence of the 
residual charge Qr on the preceding signal— the preceding sequence 
of zeroes and ones. Neither the attenuation of the size of the charge 
packet per se nor the accumulation of the incompletely transferred 
charge per se play a primary role. For example, suppose that we are 
using a very simple form of absolute-amplitude detection in which any 
charge packet of size Q > Q = (Qi + Qo)/2 is detected as a "one" 
and any packet of size Q < Q is detected as a "zero." Then a one 
preceded by a long string of zeroes will be detected as a zero Una > 0.7. 
Thus, under noiseless conditions we shall have a nonzero error rate in 
cases where the signal attenuation is only one-half (see Appendix A). 
On the other hand, the accumulation of incompletely transferred charge 
can be reduced using zero-net-charge coding. As each bit of signal is 
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coded into the same amount of total charge (distributed between two 
packets), the total amount of incompletely transferred charge from 
each bit is the same. Nonetheless, no significant improvement is ob- 
tained. Indeed, the maximum na for zero errors under noiseless con- 
ditions remains a strong function of the preceding sequence of charge 
packets (see Appendix B). This result is true for analog as well as 
digital signals, as is clear from the discussion in Appendix C. The 
dependence of the residual charge on the previous signal suggests that 
more attention should be given the detection of the signal rather than 
its coding. This we discuss in the next section. 

IV. THE OPTIMUM DETECTION OF SIMPLE DIGITAL SIGNALS 

Let us suppose that some arbitrary sequence of charge packets of 
size Qi for a "one" and Q Q for a "zero" have preceded the charge 
packet which we now wish to detect. The residue or residual charge 
added to the charge packet of interest can be designated Q R as before, 
where Q R is a function (given in Appendix A) of the preceding signal. 
If the charge packet which we are detecting is in fact a one, then the 
size Q(l) of the packet will be 

Q(l) = (1 - a )*Q x + Q R . (11) 

If, however, the charge packet is zero, then the packet's size Q(0) will 
be given by 

Q(0) = (1 - a)"Q + Q R . (12) 

One will clearly optimize the detection (even in the presence of noise) 
if one choses for the detection level Q d of the mean of Q(l) and Q(0) : 

Q d = (1 - «)"i(Q! + Qo) + Q R . (13) 

If Q > Q d , we say that we have detected a Qi packet, and if Q < Q d 
we say that we have detected a Q packet. Because Q d is a function of 
Q R which in turn depends on the entire preceding signal, we shall refer 
to this as a dynamic detection procedure. In contrast to Q d given in 
(13), the static detection procedure mentioned in Section III has 
Q> = Q/(l - a), and Q > Q, implies Qi and Q < Q, implies Q . 

It is shown in Appendix A that under noiseless conditions this scheme 
of dynamic coding is errorfree regardless of the size of na or of the 
nature of the preceding sequence of zeroes and ones. This again illus- 
trates the role of the dependence of Q R on the preceding signal, which 
we noted at the end of Section III. It is shown in Appendix A that 

0(1) - Q d = (l - «)» 2i^lQ? = Qd _ Q(0) (14) 
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as is clear from eqs. (11), (12), and (13) as well. The quantity (Qi - Q ) 
may be referred to as the dynamic range of the device. Thus relative 
to the dynamic detection level, Q d , the signal, Q(l) or Q(0), is attenu- 
ated as (1 - a) n . [Note that [Q(l) - Qi] and [Q d - Q(0)] are in- 
dependent of the residual charge Q fl .] In the presence of noise, errors 
will clearly be introduced if (1 - a) n (Qi - Qo)/2 approaches the noise 
level. This also shows that, having eliminated the signal-dependent 
residual charge, attenuation now plays an important role in limiting 
device operation. As n increases, this signal attenuation coupled with 
the compounding of noise both reduce the signal-to-noise ratio and 
lead to a reduction in the information transmission and storage 
capacities of the device. However, now it will be for not « 4 rather 
than for na tt 0.7 that attenuation becomes limiting. 

To set the dynamic detection level Q d , Qr must be realizable. In 
the absence of noise, this is always possible in principle since Q R is 
an explicit function of the known, preceding signal. In the presence 
of noise, Q R determined by eq. (31) also yields, for most cases of 
interest, nearly optimal detection in spite of the possibility that some 
preceding packets may have been incorrectly detected. 11 

In Section V, we briefly review noise in CTD's and then in Section 
VI we shall see how this dynamic detection scheme minimizes the 
error rate in the presence of noise. This further stresses the importance 
of detection in optimizing the operation of "simply coded" CTD's. 

V. REVIEW OF NOISE 

Noise in charge transfer devices is a fascinating subject which, un- 
fortunately, can be only highlighted in this section. 12-17 Owing to the 
dramatic time dependence of the current during a single charge 
transfer, the noise generated during a single transfer is quite non- 
stationary. Since nearly all theories of noise in solid-state devices 
assume that the noise is stationary, 18 it is necessary to redo much of 
the theory taking into account the nonstationary aspect. A time- 
domain analysis has been found to be most convenient, whereas 
standard treatments are carried out in the frequency domain. 

In Fig. 3 the most common sources of noise in CTD's are categorized. 
At the input, one has full shot noise only if the electrons enter the 
source independently (e.g., if generated by the random arrival of 
phonons in an imaging device or if injected by an emission-limited 
diode). At the output the nonrandom coupling to the clock line is the 
worst source of distortion in some cases. A distinction 14 is made between 
noise generated from transfer processes, typically thermal and trap- 
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Fig. 3— Sources of noise for a CTD. 

ping noise, and that from storage processes, typically leakage and again 
trapping noise. We shall return to this distinction shortly. Since a 
charge packet acquires some noise with each transfer-storage period, 
the noise component increases as the packet is transferred from input 
to output. This we refer to as compounding. Finally, the occurrence of 
clock voltage fluctuations in the presence of fabrication variations in 
the individual capacitances leads to a form of modulation noise. This 
type of noise is also compounded. For simplicity we have left out of 
Fig. 3 a number of less important noise sources. 

Let us now return to the important distinction between storage 
process and transfer process noise. 14 It should be recalled that a CTD 
shift-register performs two functions simultaneously, the transfer of 
charge and the storage of charge. In Fig. 4 we indicate the basic dis- 
tinction between the noise generated from these two processes. In the 
case of storage process (SP) noise, the charge fluctuation generated 
during each transfer cycle in each cell is essentially independent of 
that in any other cell. For transfer process (TP) noise this is not the 
case. Conservation of charge implies that if an excess of AQ is trans- 
ferred from one storage region to the next, — AQ is left behind for the 
subsequent charge packet. This introduces a correlation in the noise 
in adjacent charge packets which leads to a suppression at low fre- 
quencies of the spectral density of TP noise. SP noise, by contrast, 
is uncorrelated and, therefore, the spectral density is flat (white). 
This difference between TP and SP noise is important for analog 
applications of CTD's and is discussed in more detail elsewhere. 14 

For digital applications we shall need the ratio S/N of the square 
of the signal charge to the mean-square noise charge at the detector. 
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Fig. 4 — The distinction between storage process and transfer process noise. 

From the discussion of Section IV [see eq. (13)] the square of the 
effective signal charge at the detector is (1 - a) 2n [(Qi - Qo)/2] 2 . The 
mean-square noise charge AQ 2 can be written 

AQ 2 = AQ? nP ut(l - a) 2 " + AQiptf 8 p(n) + 2AQ 2 rP # T p(n), (15) 

where AQ 2 nput is the input noise contribution, AQ§ P the storage process 
noise acquired by a single packet during a single clock period, &Q\ P the 
transfer process noise acquired by a single packet during a single 
charge transfer, (1 - a) 2n the attenuation from input to output, 
H 8 p(n) the compounding factor 12 for storage process noise, and 
#tp(w) the compounding factor 12 for transfer process noise. These 
compounding factors are approximately equal to n for na « 1 ; how- 
ever, for na > 1 they are suppressed 12 by incomplete transfer effects. 
For na » 1, # S p(n) ^ (n/xa)* and H TF tt (2a)- 1 , both of which are 
much less than n. Essentially the effect of incomplete transfer is to 
attenuate the accumulated noise as well as the signal. Owing to the 
correlation between the transfer-process noise components, 14 #tp 
saturates. For shot noise, AQ 2 nput = qQ, where Q is the (mean) total 
signal charge (Qi - Qo). For thermal noise, AQ 2 ^ = fkTC, where T 
is the temperature of the charge carriers and C is the storage capaci- 
tance. For our purposes here we shall ignore other noise contributions. 
Thus we find 

S/N = (1 - a) 2n [(Qi - Qo)/2] 2 /AQ 2 , (16) 
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where AQ 2 is given by eq. (15). Equation (16) is plotted in Figs. 5 
and 6 for a = 10" 3 and a = lO" 4 for C = 1, 0.1, 0.01, and 0.001 pF, 
for thermal noise only and for thermal and shot noise. We shall use 
eq. (16) in Section VI on error rates and in Section VII on storage 
capacity. 

VI. MINIMUM DIGITAL ERROR RATES 

No one would operate a CTD under conditions where errors in 
detection could occur even under noiseless conditions. However, in the 
presence of noise it is possible for a "one" to acquire sufficient net 
"negative" noise charge to be detected as a "zero" even under optimum 
conditions. It is the purpose of this section to calculate the probability 
of making a detection error, and to see to what degree the error rate 
(error probability times clock frequency) is minimized for the simple, 
two-level digital coding scheme by using dynamic detection. 

Suppose that an arbitrary charge packet following an arbitrary 
sequence of charge packets would, under noiseless conditions, be of 
size Q, at the output of the shift register. In the presence of noise the 
probability P(Q) that the observed size of the packet is Q within dQ 
is given by 

P(Q)dQ = exp[-(Q - Q a y/2AQ^/(2irAQ 2 WQ. (17) 

If we are using only zeroes and ones, then the probability P of detect- 
ing a certain "one" as a zero is given by 

Pi= [ Q "P(Q)dQ, (18) 

J— 00 

where Q d is the detection level (see Fig. 7). In eq. (18), Pi depends 
upon Q, = Q(l) which in turn is a function of the sequence of signal 
charge packets preceding the one [see eqs. (11) and (12)]. To deter- 
mine the average error probability, Pi must now be averaged over all 
possible sequences of signals, in general a very difficult task. 

Let us write eq. (18) in a slightly different form by changing 
variables. 

p . fQ.*-Qa> exp(-Q72AQ*) ^ 

Pl ~ L (2.AQ*)* dQ 

or 

Pi - / e- xt ' 2 dx/(2w)* 

J— 00 

or 

Pi = /i{[Q d -Q(D]/(AQ 2 )i}, (19) 
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Fig. 5— Signal-to-noise ratio for CTD (a = 10~ 3 ) with storage capacitance C of 
1 pF, 0.1 pF, 0.01 pF, 0.001 pF. 
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Fig. 6— Signal-to-noise ratio for CTD (a = 10" 4 ) with storage capacitance C of 
1 pF, 0.1 pF, 0.01 pF, 0.001 pF. 
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Fig 7.— A charge packet of size Q. under noiseless conditions has a probability 
distribution of P(Q) in the presence of noise. The probability that Q < Qd (that an 
error is made in detection) is the area "under" P(Q) for Q < Qd- 



where, of course, 



J—a> 



(20) 



We desire (Pi) = </i>, where the brackets "< )" denote averaging 
Q(l) (and possibly Q d ) over all possible sequences of ones and zeroes 
CQ(1) = Q(l - <*) n Qi + Q R from eq. (11)]. If Q. > Qd for all possible 
Q R (as it must if errors are to be avoided under noiseless conditions), 
then we show in Appendix D that <P> ^ /{<[Q« - Q(1)]/(AQ 2 )*)}. 
In other words, the function / evaluated at the average of its argument 
is a lower bound to the average of the function, the average error 
probability. This permits putting a lower bound on the error rate. 
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If LQd — Q(l)] is independent of Q R , as is the case for the dynamic 
detection scheme discussed in Section IV, then 

<^i) = /i{[Q d - W)]/(AQ 2 )*J (21) 

and the average error probability is at the lower bound. For the 
dynamic detection scheme of Section IV, [Q d — Q(l)] = — (1 — a) n 
X (Qi — Qo)/2, which implies that 

e- xi 'Hx/{2irY , (22) 

-X 

where (S/N) is given in eq. (18). If static detection had been used 
LQd = 0/(1 — «)], then the average error probability would always 
exceed the (Pi) given in (22). 

It remains to prove that the optimum (dynamic) detection scheme 
given in Section IV gives the lowest possible error probability for 
simple digital coding. The proof is as follows. We found above that 
for detecting a one the average error probability was at least 

<i°i)ib - hlKQi) - <Q(i)>]/(AQ 2 )M- (23) 

(The "lb" stands for tower bound.) Noting the definition of /i [eq. 
(20)], we note that we can make (Pi)i b smaller by reducing (Q d ), 
(Q(l)) being already determined. However, we must also consider the 
error probability in detecting a "zero." Proceeding as for a "one," 
we obtain for P the error probability for detecting a certain "zero," 

Po= /o{[Q d -Q(0)]/(AQ 2 )M, (24) 

where now 

My) - r e-*i*dx/ '(2*)*. (25) 

J V 

Also 

<Po>i b = h\l(Q d ) - <Q(0))1/(AQ')*}. (26) 

From the definition of / , we note that we can made (P)ib smaller by 
increasing (Qd), (Q(0)) being already determined. Assuming that an 
equal number of "zeroes" and "ones" are used in the simple digital 
coding, then (by symmetry) choosing (Q d ) so that (Pi)ib = (Po)\b, 
we shall achieve the minimum lower bounds. (Pi) = (Po)ib for 
(Qd) = KQ(l)) + (Q(0))]/2. But for our dynamic detection scheme 

(Qd) = 0(1 - a)" 1 = [<Q(1)> + <Q(0)>]/2 (27) 

and for our dynamic detection scheme (Pi) = (Pi)u> and (P ) = (Po)ib- 
Therefore, since (Qd) for the dynamic detection scheme produces the 
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lowest possible lower bounds for error probabilities, and since the error 
probabilities are in fact equal to these lower bounds, no other detec- 
tion scheme can detect with lower error probability. (It is possible 
that another scheme can do just as well, however, since it is only (Q d ) 
and not Q d itself which is the determining factor.) 

It is clear that this theorem places an operational limitation (a 
minimum error rate in detection) and CTD's using simple digital two- 
level coding. The theorem can be extended 11 to the dynamic detection 
of multilevel digital codes. 

VII. MAXIMUM STORAGE CAPACITY 

One use of the CTD is as a memory or storage element. In other 
applications the CTD can be used to shift or to transfer information 
from one location to another. To properly access the operations of 
CTD's in these applications one must calculate the maximum infor- 
mation transmission capacity and the maximum information storage 
capacity of the CTD. As a result of the work of Shannon, our labors 
are greatly diminished. 

Shannon 2 proved a most profound theorem. Let B be the bandwidth 
of a transmission channel, and let S/N be the signal-power-to-noise- 
power ratio. Then the maximum transmission capacity of the channel 
Ct in bits per second is given by 

C T = B\og 2 (l +S/N). (28) 

This result can be understood for the CTD in the S/N » 1 range 
as follows. The number of levels into which a digital signal can be 
divided and still be detected with reasonably small error is (S/N)*. 
Log2(S/N) i is the maximum amount of information in bits detected 
with each charge packet. f is the rate at which charge packets are de- 
tected. Thus U log* (£/#)* = \U log 2 (S/N) W B log 2 (l + S/N) is the 
number of bits of information transmitted per second. (In Section 2.3 
we noted that for not « 1, B « f /2.) Shannon was, of course, much 
more interested in the S/N « 1 range. For this case his theorem implies 
that no matter how noisy the transmission channel may be, it is 
always possible to pass information along it. We shall not make use of 
Shannon's result in this latter range. 

A more interesting quantity from the standpoint of the CTD is the 
maximum information storage capacity. This can be calculated from 
Shannon's Theorem 2 as follows. If C T is the number of bits per second 
transmitted, then if one waits a time T equal to the time it takes the 
information to be transferred from the input to the output of the 
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linear medium, the maximum storage capacity in bits C t must be 
given by 

C. = T C T = T B log 2 (l + S/N). (29) 

For a transmission line, T is given by the length of the line divided 
by the propagation velocity. For a CTD, T = N /f where N = n/p 
and p is the number of charge transfers per clock period T = l// . 
Thus for a CTD we find for the maximum information storage capacity 
C, in bits: 

C 8 = No(B/f ) log 2 (l + S/N). (30) 

[Strictly speaking, the maximum information storage capacity will 
actually be less than or equal to the C a given in eq. (30). This is be- 
cause as S/N decreases, the length of the code word increases. 2 How- 
ever, for a CTD with N storage units, the maximum length of a code 
word is restricted to JV . Thus for small S/N, the prediction of eq. (30) 



•o" 




n NUMBER OF TRANSFERS 



Fig. 8 — Storage capacity C. versus n, the number of charge transfers, for a = 10 -3 
?( C 7 IP.volts, C = 1, 0.1, 0.01, 0.001 pF. Shown for comparison is C. for 2- and 
4-level digital codes. 
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may be in fact unrealizable. As our primary concern is for S/N » 1, 
the upper limit of C, should be sufficiently accurate.] 

Knowing N , (B/f ), and (S/N) as functions of n, the number of 
charge transfers, we can calculate C a versus n to determine the maxi- 
mum C s possible under various circumstances and for what n C„ is 
maximum. This has been done in Fig. 8 for cc = 10~ 3 and in Fig. 9 
for a = 10~ 4 . In both figures Q/C = 10 volts, Q = Qi = 2Q , and 
storage capacitance C = 1, 0.1, 0.01, 0.001 pF. Also shown is C. 
for two-level and four-level codes. Here n is limited by an n (maximum) 
for each code at the number of transfers beyond which signal degrada- 
tion due to incomplete transfer would lead to errors in absolute- 
amplitude detection in the absence of noise. We note (i) that the 
maximum C. occurs for n about a factor of three larger than for 
n (maximum) from the examples of simple coding and detection, and 
(u) that the maximum value of C, is about a factor of four to five 




n NUMBER OF TRANSFERS 



Fig. 9 — Storage capacity C, versus n, the number of charge transfers, for a = 10 *, 
Q/C = 10 volts, C = 1, 0.1, 0.01, 0.001 pF. Shown for comparison is C, for 2- and 
4-level digital codes. 
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larger than C, at n (maximum) for these common digital codes. [It 
should also be noted that C. at n (maximum) for two-level digital 
coding exceeds C, at n (maximum) for multilevel coding. This is dis- 
cussed in more detail elsewhere. 19 ] It is encouraging to note that there 
exists such a margin between what is theoretically possible and what 
can be simply accomplished. From Section VI, however, it appears that 
it will be rather difficult to achieve optimum performance if this be 
desired (or essential). 

VIII. CONCLUSIONS AND RECOMMENDATIONS 

It is not surprising to find that incomplete charge transfer and 
random noise (especially shot and thermal) limit the bandwidth, 
storage capacity, and error rate of CTD's. What is surprising, how- 
ever, is that the residual charge level Q R resulting from the portions 
of charge (preceding the packet of interest) incompletely transferred 
is so strongly signal dependent that signal detection with static detec- 
tion levels becomes seriously impaired prior to the onset of significant 
signal attenuation or noise degradation. Coding to offset the signal 
dependence of Qr is found to be ineffective for the simple examples 
considered. On the other hand, by employing our dynamic detection 
scheme, which adjusts the detection levels to null out the signal 
dependence of the incompletely transferred charge, the operational 
range is significantly extended, limited only by the physically un- 
avoidable effects of attenuation and noise. It is also shown that no 
detection scheme can be devised with a lower error rate than this 
dynamic detection scheme. 

It might be concluded on the basis of the above result that more 
attention should be focused on detection rather than coding as a 
means of offsetting the worst effects of incomplete charge transfer. 
Noting the results shown in Figs. 8 and 9, however, it is apparent that 
substantial increases in storage capacity are possible with more 
sophisticated coding-decoding schemes. 
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APPENDIX A 

In this appendix the degradation of digital signals is discussed in 
general in some detail. In Appendix B these results are applied to 
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certain specific cases using several simple coding procedures. In 
particular, the signal-dependent residual charge Qr is discussed. These 
mathematical details should be of some assistance in understanding 
several equations presented in Sections III, IV, and VI of the text. 
The analysis will be in the time domain. In Appendix C a frequency 
domain analysis is given. 

We shall denote by Q N the size of the Nth charge packet at the 
input which precedes the packet of interest by N clock cycles. As in 
the text we shall denote by Q R the size of the accumulated residual 
charge originating from the incompletely transferred portions of the 
preceding packets, Q N - Mathematically Qr is given by 7 - 19 

Q R = (l-a)"!^^)^, (3D 

where n is the total number of transfers from input to output and a 
is the coefficient of incomplete transfer for each transfer. [Equation 
(31), as well as eq. (3), are somewhat approximate. To obtain a "pure" 
binomial factor in (31), or equivalently to be able to write a single- 
transfer equation like eq. (3), one must assume that the actual transfer 
of charge can be approximated by simplified single transfers either on 
a per-cell basis as in (31) or on a storage-region basis as in eq. (3). 
The error involved in this approximation will be of the order of a or 
na 2 , whichever is larger.] 

The physical significance of (31) is the following. The portion of 
Qn which will show up in Qr are electrons incompletely transferred 
N times, each time introducing a factor a. The binomial factor gives 
the number of distinct alternative sequences of "transfer" or "no 
transfer" which can lead to a portion of Qn contributing to Qr. 

Suppose now, as in the first example in Section III, we have a 
packet of size Qi preceded by an infinite string of packets of size Q . 
Then for Qr one has 

Q«(000- ••)=(!- «)" t =i ( n ^ N ) «' v Qo 

= - (1 - a)"[l - (1 - a)- (n+1) ]Qo • (32) 

Similarly, 

g«(lll- • ) =- (1 -a)"[l - (1 - a)- (n+l) ]Qi. (33) 

For a Qi following the string of Q 's, the size of the charge packet 
Qleq. (10)] at the output will equal Q = (Qi + Qo)/2 if n is such that 

(1 - a) n Q, - (1 - a)"Q - (1 - a)-(»+ 1 >]Qo = (Qi + Qo)/2 
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or if (1 - a ) n = \ [to within a factor of (1 — a) « 1]. Should 
(1— <x) n > ?, Q < Q and the one would be detected as a zero. 

In Appendix B we shall discuss more complicated coding schemes to 
see whether Q a (31) can be reduced or at least made less sensitive to 
the signal preceding the charge packet of interest. For the present let 
us continue to derive some of the other results stated in the text. 

To cast eq. (13) into a simpler form we proceed as follows: 

Qa= (l-«)«^4^ + Q * 

= (1 - a )«Q + (1 - «)" JE (" ~$f N )a*(Qm -Q + Q) 

= [ (1 _ a)n _ a _ a)n ( 1 __L_)] +Q; 

= Q(l - a)" 1 + Q' R , (13') 

where 

Q' R m (1 - «)» J£ ( U + N y»(Q N - Q) . (34) 

The static detection level, Q, = Q/(l —a), actually differs by a 
factor of (1 - a)- 1 from the Q used in Section III and in the discussion 
following eq. (33). The difference, while insignificant, arises from 
whether one takes Q d to be the average of Qi and Q , the sizes of the 
charge packets at the input, or whether one takes Q d to be the size 
at the output of an average charge packet [of size Q = (Qi + Qo)/2] 
following a string of similar packets. Thus for such a case 

Q, = Q = Qd -«)»+ (i -"^E^^Wq 

= Q(l - a)" - (1 - o)»D - (1 - a)- (B+1) ]0 
or 

Q. = Q/(l-a). (35) 

To show that by using the dynamic level Q d given by eq. (13) one 
can have zero detection errors in the absence of noise we proceed as 
follows. Using eq. (11) one has at once that 

Q(l) -Qd' (l- a)»(Qi - Q) (36) 

independent of Qr. As Qi > Q ,Qi > Q and, therefore, Q(l) - Q d > 0. 
Similarly using eq. (12) one finds Q d - Q(0) > 0. Thus, in the absence 
of noise, Q(l) and Q(0) are always separated by Q d , and hence no error 
need be made in distinguishing them. 
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APPENDIX B 

In this appendix we use the results of Appendix A to investigate 
what improvement if any is possible in CTD operation by using several 
simple coding procedures. We shall assume noiseless absolute-ampli- 
tude detection using a static detection level at the average output 
charge level. Equation (31) for Q R can be used to calculate the result 
of other coding procedures. 

In Table I, I have enumerated four simple means of representing 
or coding a digital zero (0) and a digital one (1) using charge packets. 
The first is just to represent a by a Q packet and a 1 by a Qi packet. 
As calculated in Appendix A [see following (33)], a Qi following a 
long string of Qo's [example " (a)"] will be detected as a Q if na > 0.7. 
This is the "na Limit" entry in the table. Finally, the size of the Qi 
packet is attenuated as exp(-na) as stated. For this coding a second 
example, "(b)," is given— a • • • Q1Q0Q1Q0 • • • sequence. In this case 
Q R is always sufficiently large for a Qi and sufficiently small for a Qo 
that under noiseless conditions Q(l) > Q and Q > Q(0) for any na. 
However, as noted in Section 2.2, such a signal is attenuated as 
exp( — na) attenuation. 

One might hope that by preventing Q R from becoming much differ- 



Table I — Four Simple Means of Representing Digital 
Zeroes and Ones Using Charge Packets 



Example 


Representation 


na Limit* 


Attenuation 


(a) 1000 ■ • • 

(b) 1010 • • ■ 


1 
Qo Qi 


0.7 

oo 


e -na 
e -2na 


(a) 1000 • • • 

(b) 1010 • • • 


QoQo QiQi 


(1)0.7; (2) 1.67 
(1) 0.79 ; (2) 2.4 


e -na 
e -na 


(a) 1000 • • • 

(b) 1010 • • • 


QiQo QoQi 


(1) * ; (2) 0.8 
(1) 2.36; (2) 0.785 


e -na 


(a) 1000 • ■ • 

(b) 1010 • ■ • 


Q Qo QQi 


(1)0.4; (2)0.5 
(1)1.6; (2)3.1 


e -na 
e -na 



* The notation (1) refers to the first of the two packets forming a bit, and (2) refers 
to the second. 
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ent than Q, one could increase the na limit. In Table I three possi- 
bilities are given. The first consists merely of coding into two ad- 
jacent Q packets and 1 into two adjacent Qis. If one detects the 
second of the two Qi packets in sequence (a), the na limit is increased 
to 1.67. For sequence (b) detecting the second Qi now has an na limit 
of 2.4 while the signal is attenuated as exp( — na). How much of an 
improvement this offers, however, is questionable. To store the same 
amount of information n must be doubled reducing the 1.67 to an 
effective 0.835. To maintain the same information rate, /„, the clock 
frequency must be doubled. This will increase a : if a is doubled, then 
the 1.67 limit, already reduced to 0.835, will be reduced further to 
about 0.42. Compared with the 0.7 limit of the simplest code, this is 
rather unfavorable. One compensation is that having two packets to 
detect rather than just one can be used to reduce the error rate in- 
duced by noise. However, one can do better, as the following example 
illustrates. 

The third example in Table I is the zero-net-charge code. Here a 
is coded as a Q Q packet followed (in time) by a Qi packet. (In the 
register shifting from charge left to right this is represented as QiQ .) 
A 1 is coded as Q Qi. The advantage of this procedure is that each 
pair, whether coding a or a 1, contains the same amount of charge, 
2Q. This prevents a buildup of charge in Q R . The most demanding 
test is sequence (b) in which the na limit is 2.36. This is a significant 
improvement over the 1.67 limit in the previous example. However, 
if one takes into account that to contain the same amount of infor- 
mation n must be doubled (as now each bit requires two charge packets) 
and that the clock frequency must be doubled to maintain the same 
data rate (which will increase a), one realizes that really very little 
has been achieved by increasing the upper limit on na from na < 0.7 
to na < 2.4 (2.4/4 = 0.6). Other straightforward modifications of the 
basic 0, 1 code, of course, suffer from the same fault. Thus to achieve 
any improvement it is necessary that one still must be able to take 
advantage of the possibility of detecting both charge packets to do 
better than the simplest code. The reason for the failure of the zero- 
net-charge code in terms of frequency-domain concepts is given in 
Appendix C. 

One final example is to follow the Q or the Qi with an intermediate 
packet of size Q. As seen in Table I, sequence (a) puts an na-limit of 
0.5, which is inferior to the other codes. This attempt to reduce 
\Qb — Q\ by following Q or Q x with a Q packet to "average" out the 
incompletely transferred charge is thus seen to be ineffective. 
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APPENDIX C 

It is quite informative to briefly discuss in the frequency domain 
the effects of various digital coding schemes on the character of the 
signal. 7,9 

In Section 2.1 we noted that incomplete charge transfer leads to a 
frequency-dependent attenuation A (a>) given by 

A(a) = exp[-na(l - cos cor „)] (37) 

for /3 = in eq. (7). A(u) is plotted in Fig. 1 for various values of not. 
As discussed in Section 2.2, low-frequency components (/ « / /2) 
suffer very little attenuation, whereas components with frequency near 
half the clock frequency (ftt /„/2) are attenuated by exp(-2na), a 
large attenuation for na > 3. 

One can offset this high-frequency attenuation by the following 
scheme. If one takes every other charge packet and replaces it by a Qi 
if it originally was a Q , and by a Q if it originally was a Q h then 
relative to Q one essentially multiplies each packet in turn by +1, — 1, 
-1-1, -1, +1, -1, •••. This has the effect of converting the spectrum 
of the signal from F(f) to F(f /2 - /) : the / = component is 
attenuated as A (/ /2) and the / = / /2 component as A (0). To better 
preserve the entire signal, one can sum the outputs of a register with 
attenuation A(f) and a register with attenuation A(/ /2 - /). The 
ratio of maximum attenuation to minimum attenuation is thus im- 
proved from exp(-2na) to 2 exp(-na)/[l + exp(-2wo:)]. However, 
distortion near / = /„/4 is still significant for na > 2. 

To see the effect of the zero-net-charge coding scheme on the signal, 
consider this example. If the clock frequency is /„, the maximum fre- 
quency the CTD can carry is / /2. However, if two charge packets 
are devoted to each or 1 as in the second through fourth examples 
in Table I, then the bandwidth is reduced to / /4. If the second ex- 
ample is chosen, then the band extends from / = to / = / /4 ; if the 
third example (zero-net-charge coding) is chosen, then the band ex- 
tends from / = / /4 to / /2, the lower-frequency components of the 
signal being carried at the higher frequencies and vice versa. If 
amplified by exp(+na), the ultimate effect of incomplete transfer on a 
signal coded using zero-net-charge coding is seen to be essentially the 
same as that on a signal coded using the second example. What is most 
striking, however, is that by reducing the clock frequency by a factor 
of two and using simple coding, one reduces na by a factor of four, 
greatly reducing the attenuation. 
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By examining the effect of other coding schemes on the spectrum 
of the signal, and by taking into account the frequency-dependent 
attenuation accompanying charge transfer in CTD's it is possible to 
ascertain whether an improvement in (noiseless) detection will be in 
fact real or only apparent. 

APPENDIX D 

In noise, detection, and communication theory one often encounters 
an integral of the form 

1(A) = f A «-^/(2r)l , (38) 

J— QO 

where A > 0. This expression, while extensively tabulated numeri- 
cally, is difficult to work with analytically. In this appendix we shall 
(0 bound 1(A) between two simple analytic functions of A which 
differ by only a factor of 2, and (ii) prove that (1(A)) ^ I ((A)) for 
A ^ 0. 

(i) Bounds on 1(A) : 

In Fig. 10 we illustrate the motivation for our approximations. 
1(A) is the area under the Gaussian for x = — » to x = — A. If 
we draw a line tangent to the Gaussian at a; = — A and extend the 
line from x = — A to the x-axis as shown, the area of the triangle 
formed by this tangent, the ;r-axis, and the vertical line x = — A is 
clearly less than 1(A). Similarly, if an exponential curve [5 exp(-f-Cx)] 
also tangent to the Gaussian at x = -A and decaying to the left is 
drawn, then the area between this curve and the z-axis for x ^ - A 
is clearly greater than 1(A). Thus, if we calculate these two areas, we 
will have an upper and lower bound on 1(A). (These curves will 
clearly not cross the Gaussian if A ^ 1, the inflection point of the 
Gaussian.) 

To calculate the areas we proceed as follows. The slope of exp( — x 2 /2) 
at x = — A is A exp(-A 2 /2), and of course its value at x = — A is 
exp( — A 2 / 2). Thus the equation of the tangent is 

y(x) = exp(-A 2 /2) + A exp(-A 2 /2)(x -f A) (39) 

(which is zero for x = - A - 1/A, f(-A - l/A) = 0) and of the 
exponential is 

y(x) = exp[-A 2 /2 + A(x + A)J (40) 
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Fig. 10 — Approximating the area under a portion of Gaussian curve by bounding 
the area between that of a right triangle (whose hypotenuse is tangent to the Gaussian 
at x = — A) and by the area under an exponential (also tangent to the Gaussian at 

x A). 



Thus computing the appropriate areas we find that if A ^ 1 then 

D/2 < (2tt)*/(A) < D, (41) 



where 



D = exp(-A 2 /2)/A. 



(42) 



Such bounds are very useful in calculating error rates, where one is 
seldom interested in accuracy better than a factor of two, and where 
upper and lower bounds are often very useful. 

00 {1(A)) * I((A)), A > 0: 

To evaluate (1(A)) is clearly very difficult even under the simplest 
of probability distributions of A, whereas I ((A)) is generally very easy 
to compute if (A) is known. I ((A)) can then be used as a lower bound 
for the more interesting (1(A)). We shall now prove the above in- 
equality. (This result is reasonably well known. 20 The proof is given 
here for completeness.) 

1(A) is a function of A. According to the mean value theorem we 
may write 



1(A) = I ((A)) + ^ 



<*> 



< A -W>+5£P 



(A - (A))* , (43) 



Now then 
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where A' (A) lies between A and (A) and depends on A. Thus we may 
write 

tfU)>-I(W> + »$£L #M) (A -<*»■). (44) 

(PI 
ff2 = Aexp(-AV2)/(27r)i (45) 

which is zero or larger for A ^ 0. Thus if we are averaging A over a 
probability distribution P{A) for which P(A < 0) =0, then 
A' (A) ^ 0, and, consequently, the second term on the right-hand 
side will be zero or greater. Hence it follows that 

(1(A)) ^ I((A)). (46) 

In Section VI this inequality is used to put a lower bound on the error 
rate for detecting digital signals. 
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